Zero momentum modes in discrete light-cone quantization 
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Abstract 

We consider tlie constrained zero modes found in the apphcation of discrete hght-cone quantiza- 
tion (DLCQ) to the nonperturbative solution of quantum field theories. These modes are usually 
neglected for simplicity, but we show that their inclusion can be relatively straightforward, and, 
what is more, that they are useful for nonperturbative calculations of field-theoretic spectra. In 
particular, inclusion of zero modes improves the convergence of the numerical calculation and 
makes possible the direct calculation of vacuum expectation values, even when the zero modes are 
determined dynamically. We also comment on zero-mode contributions not included by DLCQ, 
namely zero- mode loops. 

PACS numbers: 11.10.Ef,11.30.Qc 
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I. INTRODUCTION 



The technique of discrete hght-cone quantization (DLCQ) [l|, H, H, 0] has been employed 
many times for the nonperturbative solution of various field theories, particularly in two 
dimensions. This includes recent calculations of eigenstates in supersymmetric Yang-Mills 
theories Q and (f)^ theory 0, , as well as the early applications to Yukawa theory 
(f)^ and (j)^ theories li, 9, QED (l0|, and QCD The method is based on light-front 

quantization 12| and discretization of the single-particle momentum modes. Those modes 
with zero momentum, zero modes, are not dynamical but are instead constrained by the 
spatial average of the Euler-Lagrange equation for the field 0, 13, 14 1. As such, the zero 
modes are usually neglected, because the constraint equation is considered too difficult to 
solve. 

The neglect of zero modes can have various consequences, ranging from the benign, such 
as slowed convergence of a numerical calculation, to the serious, an absence of understanding 
of vacuum effects, particularly symmetry breaking 13, 14, 15, 16, 13]. Can they instead 



be included in some straightforward way? We argue here that they can. Although the 
constraint equation can be a nonlinear operator equation, the DLCQ approximation requires 
only a finite expansion of the solution in inverse powers of the resolution K, where discrete 
longitudinal momentum fractions are measured in multiples of 1/K. The solution to the 
constraint equation can then be generated analytically order by order. For simplicity, we 
formulate the discussion in two dimensions, but the approach can be immediately extended 
to modes with zero longitudinal momentum in any number of dimensions. 

In general, the DLCQ approximation is equivalent to a numerical quadrature for the 
field-theoretic mass eigenvalue problem, where the eigenstate is expanded in Fock states with 
momentum wave functions as coefficients. The eigenvalue problem can be reduced to coupled 
integral equations for these wave functions. The quadrature points, in terms of momentum 
fractions, are equally spaced by 1/K. The zero-mode contributions are contributions from 
the numerical approximation where a momentum is zero but the integrand is nonzero, even 
though the wave function itself may be zero. To neglect such contributions is an error 
of order 1/K, which delays convergence as K oo, relative to the nominal trapezoidal 
quadrature error of l/K"^. Within the DLCQ approximation, such contributions appear 
as zero-mode contributions in the Hamiltonian. These typically take the form of effective 



interactions that include zero- mode exchange between dynamical modes |18l . |19| . Higher- 
order quadrature schemes can also generate effective interactions, which can be derived from 
the numerical approximation to the coupled integral equations for the wave functions. In 
addition, the quadrature scheme associated with DLCQ can generate contributions that 
go beyond DLCQ, to include effects from quantum corrections to the constraint equation; 
we illustrate this in the case of 0^ theory by obtaining a contribution from a zero-mode 



loop [2oi, m 



In the remainder of the paper, we develop these ideas more fully. A general discussion 
is given in Sec. [TTl followed by three specific applications in Sec. IIIII The approach is 
summarized in the final Sec. IIV[ 



II. ZERO MODES IN DLCQ 

We consider various two-dimensional theories and the contributions to their light-front 
Hamiltonians from modes of zero momentum. Much of the notation is taken from earlier 
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work on zero modes 16|. As light-front coordinates 12 
Each Lagrangian is of the general form 



we use 



(xO±xi)/v^. 



C = d< 



(2.1) 



with the scalar field of interest, V a generic interaction term which may include other fields, 
and C'^^^^'^ the free Lagrangian for any other fields. The mass term includes the possibility 
of a plus sign, in order to consider a 0*^ theory with tree-level symmetry breaking. For this 
Lagrangian, the Hamiltonian density is 



(2.2) 



We apply DLCQ [11,3 by imposing periodic boundary conditions on in the box —L/2 < 
x~ < L/2. A constant, zero- momentum mode 0o is separated from the other modes, so that 
we have (p = (p + (pQ and 



n>0 



'47m 



with = 2-nn/L and [a„, a]^] = 1. The zero mode, written 0o 



by the Euler-Lagrange equation 

(29+9_±/x2)0=_v'(0), 

which, after integration over the length of the box, yields 



=F /i 00 = 



L/2 



L/2 



V {(p + (t)Q)dx~ . 



(2.3) 



ao, is constrained [2I, 



(2.4) 



(2.5) 



This constraint is to be solved, to determine 0o and, therefore, oq in terms of the dynamical 
modes. 

The DLCQ Hamiltonian operator for evolution in light-front time x'^ is 



2L 

ax H = ±/i — 
-L/2 47r 



S2 + -Co 



1 

2' 



+ I'^^'dx- [v{<p + <p,) + n'^^r] 

-L/2 



As in 3, we define 



- y 



^n+---+in,0 



■ ai 



(2.6) 



(2.7) 



The form of E„ is made compact by using a negative index to indicate a creation operator; 
when i < 0, we have = a|^|. The expression for V~ can be simplified once the full form 
is specified and the constraint equation is invoked. Following Rozowsky and Thorn j^, we 
work with a rescaled Hamiltonian 

The eigenstates of h are constructed as Fock-state expansions at fixed light-cone momen- 
tum P+ = 271 K/L, where K is an integer that sets the resolution of the calculation [1]. The 
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momenta fc+ = 2Tm/L of the particles in each Fock state must sum to P"*", and the indices 
n must then sum to K. The hght-cone momenta, and the integer indices, must be positive; 
this hmits the number of particles to a maximum of K. 

The eigenvalue problem for h yields coupled equations for the wave functions of the Fock- 
state expansion. In this context, the zero-mode contributions come from integration end 
points where the momentum is zero. This can be tracked by starting from the continuum 
form and discretizing the integral equations for the wave functions. The discretization 
equivalent to the DLCQ approximation is a trapezoidal approximation to the integrals. In 
terms of the momentum fractions x = k^/P^, an integral from zero to one is replaced by a 
sum over discrete points in the integral, ai Xn = n/K = k^/P^, multiplied by the interval 
size, 1/K: 

f dxfix, l-x) = -^/(0, 1) + ;^ E /(^/^' 1 - ^/^) + ^/(l' 0) + V^')- (2.9) 

"^0 n=l 

The end-point corrections, which are the zero- mode contributions, are then of order 1/K 
higher than the bulk of the sum that approximates the integral. Thus, the zero-mode 
contributions are equivalent to the addition of effective interactions to the Hamiltonian, 
interactions that include explicit powers of 1/K [3], and have the effect, at a minimum, of 
improving numerical convergence by restoring the 1/K terms missed in ignoring the end- 
point corrections of integrals. 

In DLCQ these terms are generated by the zero-mode part of the field as determined 
by the constraint equation. For consistency with the underlying trapezoidal approximation, 
these contributions should be kept to no higher in 1/K than the corrections expected. The 
order is measured relative to the non-zero-mode parts, which are typically proportional to 
S„ = where 

Sn = ^ E 4^=^:a,---a.„: (2.10) 

is written explicitly in terms of momentum fractions Xj. A non- zero- mode contribution of S„ 
would then require zero-mode contributions of no more than order l/i^^+"/^. We illustrate 
this in the following section for various interaction models. 



III. APPLICATIONS 

A. Wick Cutkosky Model 

The simplest nontrivial case is an interaction V = X(f)\x\'^ with a complex scalar field x of 
mass m. The spectrum of this theory is unbounded from below [22|; this is obvious classically 
because the field can acquire a negative value and an arbitrarily strong x field can then 
drive V to —oo. An ordinary DLCQ calculation that excludes zero modes can detect this 
structure, but not without careful extrapolation [23|]. Here we include zero modes, and a 
simple variational calculation is all that is needed to detect the unbounded behavior. 

We use antiperiodic boundary conditions for the discretized x field and thereby avoid its 
zero modes. This option is not available for the field, since it is coupled to the necessarily 
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periodic square of the x field. The mode expansion of the x field is 

X = Y^ ke^'^"^" + 4e-*^"^"j . (3.1) 



^„ V L 

n>0 * 



The constraint equation (12. 5p yields 

ao = -2^S^, (3.2) 
^2 = ^ E - K^n + 4^^n] . (3.3) 



where (? = A//i v47r and 



The Hamiltonian reduces to 

h = \^2 + - / {T^lf + %==«^ + ddm] , (3.4) 

where the sum in the last term includes both positive and negative indices. The (Sg)^ 
term has as its origin the zero- mode contributions (pQ and 0o|xP to the Hamiltonian density. 
The constraint equation specifies that (po is proportional to the negative of the average of 
IxP, making a net negative-definite contribution to the energy density. 

To isolate this negative contribution, consider the expectation value of h for the highest 

Fock state ^c|j (^d\^ |0) of K x-particles in any charge sector, each with the same 
momentum fraction 1/K. The expectation value is 



Since this tends to —00 as K ^ 00, the spectrum extends to —00 in the continuum limit. ^ 



B. 0^ Theory 

In 0^ theory there are, of course, no fields other than 0, and the interaction is = ^0^. 
The constraint equation (12.51) becomes 

- ao = + (3.6) 

with g = A//i^-\/47r. On use of this constraint, the Hamiltonian can be written 

h = \^2^ - ^al - ^al. (3.7) 



This fate is avoided in Yukawa theory, where x is a fermi field, simply because the identical x — \/K 
states cannot be populated. 
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The non-zero-mode piece of the Hamiltonian, the first two terms, can be written as ^2 + 
^^^3/2 S3, where the S„ are defined in Eq. (12.101) . Thus the zero-mode contributions we seek 
are of order l/K^^"^ at most. 

The constraint equation is solved to a consistent order by taking = vq + Vi/K and 
finding vq and vi. From the constraint equation we have 

- ^0 = l^^o a^nd - f 1 = VqVi + ViVq + gT,2- (3.8) 

For vq, the two possible solutions are Vq = 0, the local minimum in the classical energy 
density, and Vq = —2/g, the local maximum. Either solution is acceptable as a starting 
point, but the first leads to simpler expressions and is sufficient for our purposes. The 
second corresponds to a different choice for the perturbative vacuum, annihilated by the 
dynamical a„, because it represents a different choice for (p = (p — (pQ] however, this change 
in the perturbative vacuum is compensated by different terms in the Hamiltonian, with no 
net effect on the spectrum. With vq = 0, the solution for vi is then immediately Vi = — (7S2, 
and we find for the zero mode 

ao = -|s2 + 0(^). (3.9) 
On substitution of this expansion for the zero mode, the Hamiltonian becomes 

The expectation value of h for the state populated with K particles of momentum fraction 
1/K is 

W = f-f^^'- (3.11) 

This tends to —00 in the K —>■ 00 limit, and, as known for a cubic theory j22|, the spectrum 
is unbounded from below. 



C. 4)^ Theory 

For 0^ theory, with its interaction V = ^(p'^, the contributions of zero modes are more 
subtle than for the first two applications. Their inclusion will improve numerical convergence 
and may provide a means to understand vacuum structure and symmetry breaking. Previous 
calculations [g], 0] of the spectrum did not include zero modes, and previous studies of the 
constraint equation [l^ attempted to solve it fully, rather than keeping ao only to an order 
consistent with the DLCQ Hamiltonian. 

The constraint equation (12.51) in this case is 

=F ao = f ao + 2gT,s + ^ (aoSs + Ssao) + f -r^a„aoa_n, (3.12) 
3 3 3 ^-^ \n\ 

where g = X/Sirfi"^. The Hamiltonian is 

J_ 



n - ±^^2 + -^^4 - — ao + 24^3/2 2^ , W^QoQm + afeaoa^a. 



cSl?y^.~~^^riala~n- aQCLna^naii)- (3.13) 

Z4A ^—^ X\n\ 
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From the K dependence of the non-zero- mode pieces of the Hamiltonian, we see that we 
need zero-mode corrections to order and thus expand cq to order l/K^^"^. 

To solve the constraint equation to this order, we write oq = vq-\-vi/ K^f'^+V2/ K+v^/ K'^^'^ 
and solve for the Uj. Two of the three possible solutions are available only for the wrong-sign 
mass term. They have Vq = ^^/sjg and correspond to the local minima in the classical 
energy density. The third solution, vq = 0, yields the local maximum as well as simpler 
expressions; it is also the only solution in the case of the correct-sign mass term. For 
simplicity, we choose to work with this third solution, for which vi and V2 are also zero and 

«o = T-^Ss = t29^3. (3.14) 

Here the upper (lower) sign corresponds to the upper (lower) sign in the mass term of the 
Hamiltonian. 

In the Hamiltonian, this solution to the constraint equation, combined with the keeping 
of terms to consistent order, yields 



^2^^' + "!^^^^ tIr^ 5Z -^==i^kaif:3am + ak^saia^) + 0{l/K^). (3.15) 



The first term is the (rescaled) mass term, and the second is the ordinary interaction term. 
The third term provides the zero-mode corrections to the mass eigenvalue problem at an 
order in 1/K that is consistent with the DLCQ approximation. The sign of the term is 
particularly significant. In the case of a positive mass term, this correction is negative, 
which will allow the spectrum to extend below zero and provide a nontrivial vacuum state 
and symmetry breaking at larger couplings. For the case of a negative mass term, where 
the symmetry breaking effect is built in, the zero-mode term is positive and can play a role 
in restoring the broken symmetry and the trivial vacuum at larger couplings. 

It is known 0, iH that for (j)'^ theory there is an additional 1/K correction from a 



zero-mode loop. Such loops do not enter into DLCQ but were found in discretization of 
perturbative contributions to forward scattering [lO]. They can also be found from a non- 
perturbative perspective by considering the coupled equations for the n-particle Fock-state 
wave function ipn{xi, . . . , Xn), which are 

1 ^ _|_ ^7/3 'lpn-2{Xl, ■ ■ ■ , Xi + Xj + Xk, ■ ■ ■ , Xn) (3 16) 

^ " ^yn{n - 1) ^JxiXjXk{xi + Xj + x^) 



dX-^^d,X2'lpn+2i-^l: • • • ; -^l? • • • j ■^2' • ■ • i -^i -^1 -^2' • • • ; •^n) 



g \ f dx Ipni^l: ■ ■ ■ , X , . . . , Xi -\- Xj X , . . . , Xn) / n r2 / 2\ i 

+ 2 / / Ht = 

J Y XiXjX yXi 'T Xj X ) 

One of the contributions to ipn+2{.yi) is 

5'/3 1 ipniyi,---,yi + yj + yk,---,yn) 

MV/i' - Er""' j- v/(^ + 2)(n+l) ^y^y^y,^y^ + y^+y,) ' 
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and substitution into Eq. fl3.16p yields a term of the form 



IE 



XiX'2Xi[Xi Xi X'2) ^^-^ / L^j^i x\ xL 



V^n(2'l) ■ ■ ■ i^i 3^1 2^2; • • • ) ~l~ -^l ~l~ •^li • • • 1 



■^x-|^X2a;fc(x-|^ + X2 + Xk) 



which represents the exchange of two bosons, with momentum fractions x[ and Xg, between 
the ith and fcth particles of the n-body Fock state. In the corner of the two-dimensional 
integral, where both x'^ approach zero and therefore two zero modes are exchanged, the 
integrand is singular. However, this is an integrable singularity; the net effect is for there to 
be a contribution of order 1/K from this corner, instead of the nominal order of l/K"^, which 
DLCQ would neglect. Where only one approaches zero, which represents an exchange 
of one zero mode and one non-zero mode, we have a 1/K contribution that is included by 
DLCQ. This contribution is generated by the new third term in the Hamiltonian fl3.15p . A 
term to represent the zero- mode loop contribution can also be added to the Hamiltonian [l^] . 

In addition to the zero-mode corrections to the Hamiltonian, the zero mode also allows 
for a direct calculation of the DLCQ approximation for the vacuum expectation value. The 
eigenstates of h separate into sectors of odd and even particle number, and in the continuum 
limit the lowest state of each sector become degenerate. Let |o) and |e) be the odd and even 
ground states and form maximally mixed states |±) = (|o) ± |e))/ ^/2. The expectation value 
for the field is then 

... ^ ^ . ±». (3.19) 

Since changes particle number by one, only the cross terms contribute to the numerator, 
and since the states have the same momentum, only the zero mode of the field can contribute. 
Thus, we have 



(e|0o| 


|o) + (o 


I0o| 


\e) 


(e| 


|e) + (o| 


\o) 



(3.20) 

On substitution of the solution (13.14p to the constraint equation, this becomes 

(3.21) 



^ 9 (e| 


IS3I 


|o) + {o\ 


IS3I 


\e) 




(e 


|e) + {o\ 


\o) 



So, one can calculate v if the mass eigenvalue problem is solved to find the lowest odd and 
even states, and one can study the continuum limit as the resolution K ^ 00. 



IV. SUMMARY 

Although zero modes are traditionally neglected in DLCQ calculations, they can actu- 
ally be taken into account without much additional effort. Each mode satisfies a constraint 
equation that is the (light-front) spatial average of an Euler-Lagrange equation. The con- 
straint connects the zero mode to the dynamical degrees of freedom and can be solved either 
explicitly or in terms of an expansion in inverse powers of the DLCQ resolution K. The 
latter expansion is truncated at the order consistent with the DLCQ approximation. 
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From the zero-mode solution, one can construct contributions to the Hamiltonian which 
will, at a minimum, repair the convergence of the DLCQ calculation. The DLCQ approx- 
imation is essentially equivalent to a trapezoidal approximation to integral equations for 
Fock-state wave functions, with quadrature points spaced equally by 1/K. The error for 
such an integral is of order XjK"^^ but when the zero modes are neglected the error Y&\jK 
and convergence as ii" oo is slowed. The zero-mode contributions to the Hamiltonian 
repair this and restore the behavior of the integration errors. For some theories, 

such as 0"^ theory, there are additional corrections beyond DLCQ but still of order iZ-fT, 
such as zero-mode loops, that correspond to order-/i quantum corrections to the zero-mode 



Inclusion of zero modes also makes possible the direct calculation of vacuum expectation 
values. Explicit symmetry breaking will yield a c-number contribution to the zero mode 
which trivially has a vacuum expectation value, but this is not the only way to have a 
nonzero value. The zero mode will, in general, have dynamical contributions. For these, 
the expectation value will be zero with respect to the perturbative vacuum. However, the 
spectrum of the theory can be such that the perturbative vacuum is not the eigenstate of 
lowest energy. Instead, some nontrivial eigenstate has a lower energy and, with respect to 
this state, the dynamical contributions to the zero mode can have a nonzero expectation 
value. 
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